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Abstract 

A measurement of charm and beauty photoproduction at the electron proton collider HERA 
is presented based on the simultaneous detection of a D*^ meson and a muon. The cor¬ 
relation between the D* meson and the muon serves to separate the charm and beauty 
contributions and the analysis provides comparable sensitivity to both. The total and dif¬ 
ferential experimental cross sections are compared to LO and NLO QCD calculations. The 
measured charm cross section is in good agreement with QCD predictions including higher 
order effects while the beauty cross section is higher. 
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1 Introduction 


At the electron proton collider HERA, heavy quarks are predominantly produced via photon- 
gluon fusion, 'jg cc or bb, where the photon is emitted from the incoming lepton and the 
gluon from the proton. The production cross sections are largest for photoproduction, i.e. for 
photons with virtuality — 0. The light quarks u, d and s are produced much more copiously 
than c and b, and beauty production is suppressed by a factor of approximately 200 compared 
to charm. Charm and beauty measurements performed at HERA so far relied on the tagging of 
only one heavy quark in each event. While the charm measurements [1-11] were mostly based 
on the reconstruction of D mesons, the beauty measurements [12-17] used semi-leptonic decays 
or lifetime signatures or both. Here an analysis is presented, where in a large fraction of events 
both heavy quarks are tagged using a D*^ meson and a muon as signatures. The correlations 
between the direction of the muon with respect to the D*^ and their electric charges are used 
to separate charm and beauty contributions. Total cross sections are measured separately for 
the processes ep eccX —> eD*pX' and ep —ebbX eD*pX' in the visible kinematic 
region, while differential cross sections are derived for combined samples of cc and bb events. 
The measurements, which are based on an integrated luminosity of C = 89 pb“^, are compared 
to leading order (LO) and next-to-leading order (NLO) perturbative QCD (pQCD) calculations. 

This measurement extends to significantly lower centre-of-mass energies of the bb system 
than previous measurements of beauty cross sections at HERA. The simultaneous detection of 
the D* meson and the muon makes possible new tests of higher order QCD effects. Eor instance, 
in the photon-gluon rest frame the angle between the heavy quarks is 180° at leading order, but 
at next-to-leading order it can differ significantly from this value due to hard gluon radiation. 
Eurthermore, the D*p pair is expected to be sensitive to a possible transverse momentum kt of 
the gluons entering the quark pair production process. 


2 Separation of Charm and Beauty 

The separation of charm and beauty contributions exploits the charge and azimuthal angle' 
correlations of the D* meson and the muon. The azimuthal angle difference A® between the 
D* and the muon and their respective electric charges Q{D*) and Q{p) are used to define four 
‘correlation regions’ I-IV. Eor Q{D*) = Q(/i) regions I and II cover A® < 90° and A<l> > 90°, 
respectively. Regions III and IV are defined correspondingly for Q{D*) ^ Q{p). 

The four regions are populated differently by charm and beauty events as is illustrated in 
figure [H Neglecting any transverse momenta of the photon and the gluon, the fusion process 
75 f —cc or bb leads to a back-to-back configuration of the two heavy quarks. Approximating the 
directions of the D*^ meson and the muon with those of the quark and antiquark and selecting 
opposite charges, cc pairs populate correlation region IV. In contrast, beauty events populate 
regions II, III and IV, depending on whether the muon originates from the same b quark as the 
D* or from the opposite b. If the muon originates from the same b quark as the D* meson, 

'The coordinate system is defined with the z-axis pointing in the proton beam direction and x (y) pointing in 
the horizontal (vertical) direction. The azimuthal angle $ is measured in the x-y plane and the polar angle 9 with 
respect to the z direction. 
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the events lie in region III. For muons eoming from the b opposite to the D* meson, the direet 
deeay populates region II, while the easeade proeess h ^ c ^ ii populates region IV. Region 
IV hence receives contributions from both cc and bb events and region I stays empty. 

The azimuthal angle correlations are smeared by fragmentation and semileptonic decay pro¬ 
cesses and by higher order QCD effects such as gluon radiation and any initial transverse mo¬ 
mentum of the gluon. Processes such as heavy quark decays via r leptons conserve the charge 
correlation, which is not the case for B^-B^ mixing and e.g. the decays b —> cW~; W~ —*• cs. 
According to the PYTHIA Monte Carlo simulation [18], which takes into account these smear¬ 
ing effects, the relative population of the four regions is as given by the numbers in figure d 
These numbers apply in the analysed kinematic region defined in section IdAl and do not include 
detector effects. Since the population of the four correlation regions is very different for bb and 
cc, it can be used for the separation of these two components. 



Figure 1: Definition of the correlation regions in terms o/A$ and the relative charges of the 
D* meson and the muon. The sketches illustrate these correlations in cc and bb quark decays 
to D*ti. The numbers represent the relative distribution over the four correlation regions of 
charm and beauty events that satisfy the cuts given in section WT\ as obtained from the PYTHIA 
simulation without detector effects. 
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3 QCD Calculations and Monte Carlo Simulations 


The Monte Carlo simulations PYTHIA [18] and CASCADE [19] are used for the deseription 
of the signal and background distributions in the separation of charm and beauty, for the de¬ 
termination of efficiencies and acceptances and for systematic studies. Their predictions are 
also compared with the measured cross sections. In PYTHIA and CASCADE leading order 
matrix elements which take into account the mass of the heavy quark are implemented and 
parton showers in the initial and final state are included to approximate higher orders (EO- 
ME-i-PS). The parton evolution in PYTHIA uses the DGEAP equations [20]. In addition to the 
direct process, a resolved photon component is generated in PYTHIA where the photon fluc¬ 
tuates into a hadronic state acting as a source of partons, one of which participates in the hard 
interaction. This component is dominated by heavy flavour excitation processes [21], where 
the heavy quark is a constituent of the resolved photon. In the PYTHIA calculation of heavy 
flavour excitation, in which quark masses are neglected, the contribution of excitation to the 
total charm cross section in the analysed kinematic region is found to be 41%, while it is 23% 
for beauty. In comparison to this component the contribution of the resolved component due to 
light quarks or gluons in the photon can be neglected in the present analysis, as can the heavy 
flavour component of the proton. 

CASCADE contains an implementation of the CCEM [22] evolution equation for the ini¬ 
tial state parton shower. The 'yg —cc or bb is implemented using off-shell matrix elements 
convoluted with kt unintegrated proton parton distributions. PYTHIA and CASCADE use the 
JETSET program as implemented in PYTHIA for the hadronisation (via the Peterson fragmen¬ 
tation function [23]) and for the decay of beauty and charm quarks. In order to correct for 
detector effects, the generated events are passed through a detailed simulation of the detector 
response based on the CEANT program [24] and the same reconstruction software as used for 
the data. 

The measured cross sections are also compared with NEC pQCD calculations in the massive 
scheme [25] using the program EMNR [26]. These calculations are expected to give reliable 
results in the kinematic region considered here, where the transverse momentum of the heavy 
quark is of the same order of magnitude as its mass. The calculations are available for both 
the direct and resolved photon processes. However, in contrast to the PYTHIA program, heavy 
flavour excitation is not explicitly included in the resolved part of the EMNR program. The 
contributions of the resolved light quark and gluon components are found to be small in EMNR 
(< 3% for charm and < 6% for beauty in the analysed kinematic region) and are neglected. 

The original EMNR program is extended to include the effects of the hadronisation of c and 
b quarks and their semileptonic decays in order to make comparisons with the measured cross 
sections in the experimentally accessible kinematical region. The heavy quark is ‘hadronised’ 
by rescaling the three momentum of the quark using the Peterson fragmentation function. Eor 
the decay into muons the momentum spectrum is implemented as obtained from JETSET. In 
the case of beauty quarks, the direct decays of 6-fiavoured hadrons into muons are taken into 
account as are the decays via a charm quark, 6 —> c —> /i. When the D* meson and the muon 
originate from the same quark, the angular and momentum correlations are implemented as in 
JETSET. The measured fragmentation fractions [27,28] for c and b quarks given in table [T] are 
used for the calculation of the cross sections. The important parameters of the pQCD programs 
used in this analysis are summarised in table|2l 
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Fragmentation Fractions [%] 


23.5 ±0.7 
17.3 ± 1.6 
10.03 ±0.64 


c 

b 

b 


/i 

±>*/i 


9.8 ±0.5 
10.95 ±0.27 
2.75 ±0.19 


Table 1: Fragmentation fractions [27, 28] used in the QCD calculations of the cross sections. 
The b ^ c ^ fi fraction also contains the b ccs decay and the r contributions. 



PYTHIA CASCADE FMNR (NLO) FMNR (LO) 

Version 

6.1 1.2007 

Proton PDF 

Photon PDF 

CTEQ5L [30] J2003 [19] CTEQ5M [30] CTEQ5L [30] 

GRV-GLO[31] - GRV-GHO[31] GRV-GLO[31] 

Renorm, scale fi] 
Factor, scale qf 

+ PTq 4m2 + ml + 

ml +PTq s + Qt ml + pl,gg, A{ml + p|.^g) 

mb [GeV] 
me [GeV] 

4.8 4.8 4.75 

1.5 1.5 1.5 

Peterson eb 

Peterson ec 

0.008 0.008 0.0033 (0.42 for 679* ) 0.0069 

0.078 0.078 0.035 0.058 


Table 2: Parameters used in the Monte Carlo and NLO programs. The FMNR calculations are 
performed in the MS scheme using the default values of Kqcd for the parton density functions. 
Pr ond pf denote the renormalisation and factorisation scales, p^^^ the average of the squares 
of the transverse momenta of the two heavy quarks, rriq the heavy quark masses, s the centre- 
of-mass energy squared, the transverse momentum squared of the heavy quark system, pxq 
the transverse momentum of a heavy quark and Cg the Peterson fragmentation parameters [29]. 
For the B^-B^ mixing values of Xd = Am^o /F^o = 0.73 and Xg = 18/257 are used. 


4 Data Analysis 


The data were collected with the HI detector [32,33] at HERA during the years 1997 to 2000 
and correspond to an integrated luminosity of £ = 89 pb”\ The largest part of the luminosity 
(80%) was collected at a centre-of-mass energy of ^/s ^ 320 GeV, the beam energies being 
27.6 GeV and 920 GeV for electrons^ and protons, respectively. The remaining 20% of the 
luminosity was taken at ^/s ^ 300 GeV (proton energy 820 GeV). 


4.1 Event Selection 

A detailed account of this analysis can be found in [34]. Events with at least one reconstructed 
D* and at least one muon are selected; multiple D* or muon combinations are treated as separate 
events. The D* is reconstructed via the decay channel^ D*~^ —> D^irf —> K~7r'^Trf (branching 
ratio (2.59 ± 0.06)% [28]), where Tig refers to the low momentum vr in the decay. The decay 

^HERA has been operated with electron and positron beams. These periods will not be distinguished in this 
analysis. 

^Charge conjugate states are always implicitly included. 
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particles of the D* meson are reeonstrueted in the eentral traeking deteetor (20° < 6 < 160°) 
without partiele identifieation. Muons are identified by reeonstrueting traek segments in the 
instrumented iron return yoke of the solenoidal magnet. These are linked to tracks in the eentral 
traeking deteetor. In order to ensure good deteetor aeeeptanee, outs on the transverse momentum 
Pt with respeot to the proton direetion and the pseudorapidity = — In tan(6*/2) are applied 
for the D* meson and the muon in the laboratory frame (see table QJ. 

Photoproduotion events are seleeted by demanding the absenee of any signals for the seat- 
tered eleetron, restrieting the aooepted range of negative four-momentum transfer squared to 
be below 1 GeV^. A out on the inelastioity 0.05 <y<0 .75 is applied, where y = P ■ q/P ■ k 
(q, k and P are the four veotors of the exehanged photon, ineoming eleetron and proton, re- 
speetively). The variable y is reeonstrueted from the measured hadronie final state using the 
Jaoquet-Blondel method [35]. The events are triggered by fast signals from the eentral traeking 
and muon deteotors. The analysed ‘visible’ kinematio region of the measurement is defined in 
table |3l 


Seleetion of D* D^tTs —> KtitTs 

Pt{K),pt{tt) > 0.4 GeV 

Pt{71s) > 0.12 GeV 
mK-jT — rnoo \< 0.080 GeV 

AM = mK-KTTs ~ itT'Ktt < 0.1685 GeV 

Visible kinematie region 

Pt{D*) > 1.5 GeV 

1 p{D*) \< 1.5 
p{p) > 2 GeV 

1 p{p) \< 1.735 

0.05 <y < 0.75 

Q2 < 1 GeV^ 


Table 3: The D*selection cuts and definition of the visible kinematic region. 


4.2 Fit Procedure 

The D* yield is measured using the AM teohnique [36], where AM = mK-wKs ~ 'i^Kn is 
the difference of the invariant masses of the KtitIs and the Ktt systems. Figure |2^ shows 
the AM distribution for the seleeted D*jji sample separately for the ‘right’ (JT“7r+7r+) and 
‘wrong’ {K~Ti~Tif) oharge eombinations. The wrong eharge distribution is normalised to the 
right eharge distribution in the range 0.155 < AM < 0.1685 GeV. The number of signal 
events is extraeted from a fit to the AM distribution using a Gaussian funetion for the signal 
and a parameterisation of the baekground"^. The parameters of the baekground funetion are 
determined from right and normalised wrong eharge eombinations. The result of this fit for the 
total signal is also shown in figure 

The total number of D*y events obtained from the fitisA^ov = 151±22. This number 
still eontains a eontribution from ‘fake muons’, i.e. from hadrons misidentified as muons and 

■^The functional form used is ci (AM — rriTrY^ (1 ~ ^3 (AMf) where Ci are fit parameters. 
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muons from the decay of light mesons. The fake muon background contributes about 37% 
in charm and 5% in beauty initiated events according to the respective PYTHIA Monte Carlo 
simulations, which give an adequate description of the data. 

In figure |2b-e> the AM distributions of the selected D*n events are shown separately for 
the four correlation regions defined in section |21 Clear peaks due to 77* mesons are observed in 
regions ITIV, whereas region I shows little or no signal, consistent with the expectation. 

The charm and beauty contributions in the data are determined by performing a simulta¬ 
neous likelihood fit of the AM distributions in the four correlation regions. In the following, 
this fit will be referred to as a ‘two-dimensional fit’, in order to distinguish the results from the 
separate one-dimensional fits of AM in each correlation region. 

In this two-dimensional fit, in addition to the D*fi contribution from bb and cc, the fake muon 
background and the combinatorial background under the AM peaks have to be considered. The 
position and width of the AM peak corresponding to the D* signal as well as the parameters 
describing the shape of the combinatorial background are fixed to the values obtained from the 
one-dimensional AM fit to the total sample (figure|2K)- The normalisation of the combinatorial 
background is fitted using right and wrong charge combinations in each region separately. The 
relative distributions of signal events from charm and beauty between the correlation regions as 
well as the fractions of fake muon background in each region predicted by the PYTHIA Monte 
Carlo simulations are used as input for the fit. In total there are six free fit parameters, the total 
numbers of D*n events from bb and cc quark pairs, Nf, and Nc, and four parameters for the 
combinatorial background, one in each correlation region. 


5 Results 


The result of the two-dimensional fit is shown together with the data in figures |2j5-e. The data 
are described well and the quality of the fit is good (x^ = 145.5 for 154 d.o.f.). In figure |3l 
the numbers of D* signal events from the two-dimensional fit in the four correlation regions 
are compared to the results of one-dimensional fits of the AM distributions performed in each 
correlation region separately. The agreement is very good. The distribution of the contributing 
processes as obtained from the two-dimensional fit is also shown in figure |3l The following 
event numbers and errors are obtained for the charm and beauty contributions from the two- 
dimensional fit: 


Ac = 53 ± 12 Aft = 66 ± 17. 

5.1 Total Cross Section 

The number of b and c events are used to compute the total cross sections in the ki nematic 
region defined in table Q] The efficiencies and acceptances are derived from the Monte Carlo 
simulations. Values of 
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Charm 

Cross section [pb] 

Data/Theory 

Data 

250 ± 57 ± 40 


PYTHIA (direct) 

242 (142) 

1.0 

CASCADE 

253 

1.0 

EMNR 

286^1^ 

0.9 

Beauty 



Data 

206 ± 53 ± 35 


PYTHIA (direct) 

57 (44) 

3.6 

CASCADE 

56 

3.7 

EMNR 

521^^ 

4.0 


Table 4: Measured D*^ cross sections for charm and beauty production in the kinematic re¬ 
gion defined in table E] For the data the statistical and the systematic errors are given. The 
LO-ME-vPS predictions (PYTHIA, CASCADE) and NLO calculations (FMNR) are also shown. 
The uncertainties of the FMNR results are obtained by varying the renormalisation and the 
factorisation scales simultaneously by factors of 0.5 and 2. The uncertainty due to a variation 
of the quark masses nic by ±0.2 GeV and nib by ±0.25 GeV is added quadratically. The last 
column shows the ratios of the measurement to the prediction. 


nl^fiep —>• eD*p,X) = 250 ± 57 (stat.) ± 40 (syst.) pb 

and of 

'^visi^P e-D*p,X) = 206 ± 53 (stat.) ± 35 (syst.) pb 

are obtained for charm and beauty production, respectively. The measured cross sections are 
similar due to the definition of the visible kinematic region, which requires in particular a high 
momentum muon, suppressing central charm production. The results are compared with the 
pQCD predictions in table S where error estimates due to the uncertainty of the quark masses 
and the scales are given for the NLO calculations. In order to assess the influence of mass 
effects in the extraction of gluon densities used in the calculations, the default CTEQ5M sets 
have been replaced by the CTEQ5F sets [30]. The results are found to be compatible. The 
uncertainties for PYTHIA and CASCADE are found to be of similar size as those of the NEO 
calculations. The measured cross section for charm production agrees well with the EO-ME+PS 
models (PYTHIA and CASCADE) and the NEO prediction (EMNR). The measured beauty 
cross section exceeds the calculated cross sections. In other recent measurements of the beauty 
cross section in photoproduction at HERA [12,15], based on the selection of high transverse 
momentum jets, ratios of measurement and EMNR based calculations between 1 and 3 are 
found. Note that the present analysis extends down to the production threshold for bb, while the 
jet measurements have a threshold which is approximately 5 GeV higher in the bb centre-of- 
mass system. 

The systematic uncertainties of the cross section measurement are evaluated by varying the 
Monte Carlo simulations. The dominant experimental errors come from the uncertainties in the 
track reconstruction efficiency (13%), the trigger efficiency (5%) and the width of the AM sig¬ 
nal (3%). Smaller contributions are due to uncertainties in the determination of the background 
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due to misidentified muons^ (1%;1.5%) and in the fragmentation fractions (1%;1.5%). Model 
uncertainties are estimated using the CASCADE Monte Carlo generator instead of PYTHIA 
(3.5%) and either taking into account or omitting the resolved component in PYTHIA (3%;5%). 
Taking into account the uncertainties due to the contribution of D* reflections (5%), the muon 
identification, the luminosity measurement and the D* decay branching ratios, the total system¬ 
atic errors for the charm and beauty cross sections are estimated to be 16% and 17%, respec¬ 
tively. 


5.2 Differential Cross Sections for Charm and Beauty 

Differential cross sections for D*n production in the visible kinematic region are evaluated as 
functions of variables characterising the D* meson, the muon and the D*n system. In this 
section results are presented for the complete data set, which contains the contributions from 
charm and beauty (figures El and |5l). 

In order to compute the differential cross sections for the data, the numbers of events in bins 
of the chosen variable are determined by a fit to the AM distribution in each bin, as described in 
section lO Here, no attempt is made to separate charm and beauty contributions. A correction 
for ‘fake muons’ is applied according to the Monte Carlo simulation. Since the fake muon 
fraction is different for charm and beauty, it is computed using the b fraction of 45% given by 
the measured cross sections (table El- 

The data are shown with the results of the PYTHIA and CASCADE Monte Carlo models 
and the EO and NEO EMNR calculations. In the theoretical models, the beauty and charm 
contributions are combined according to the measured total visible cross sections (table El and 
normalised to the sum of these cross sections, in order to facilitate a shape comparison. The 
error bands for the EMNR prediction are computed as for the total cross section (see caption 
of table El). The measured differential cross sections are similarly normalised, which has the 
advantage that the systematic errors largely cancel and are negligible compared to the statistical 
errors. 

EigureElshows the differential cross sections as a function of the transverse momentum and 
pseudorapidity of the D* meson and the muon separately. Overall the QCD models describe 
the shapes of the measured distributions quite well, although there is a tendency for the mea¬ 
sured pt{D*) and Pt{iA distributions to be softer than the calculations. A slight discrepancy 
is also present in the differential cross section as a function of the pseudorapidity of the muon 
(figure Efc) which shows a central dip due to the large muon momentum required. 

Quantities derived from a combined measurement of the D* and muon are shown in figure|5l 
In figures |5t and c, the differential cross sections as a function of pt{D*p), which is defined as 
Pt{D*p) = \pt{D*) +Pt{p)\, and A$ are compared with the EO and NEO EMNR predictions. 
The data show the expected deviations from the EO calculations due to higher order effects: the 
observed pt{D*p) distribution is flatter and the A$ peak around 180° is broader than the EO 
computation. The data are in good agreement with the NEO calculation. In figures |5j5 and d, the 
same differential cross sections for pt{D*p) and A® are compared with PYTHIA and CAS¬ 
CADE which also give a good description of the data. Although different approaches are used 

^Where two numbers are given the first applies to charm and the second to beauty. 
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in PYTHIA and CASCADE to compute the evolution of the partons from the proton and the 
hard interaction, the differences between the two simulations are smaller than the experimental 
errors. 

Figures |5t and f show the invariant mass, M{D*n), and the rapidity^, y{D*n), of the D* 
meson and the muon together with NLO FMNR, PYTHIA and CASCADE predictions. The 
invariant mass M{D*jj) reflects the eentre-of-mass energy of the quark pair and y{D*jji) is 
related to the ratio of the energies of the partons entering the hard interaetion from the proton and 
the photon. Both differential cross sections are adequately described by all model calculations. 


5.3 Results for a Charm Dominated ‘Quark Antiquark Tag’ Sample 

The cross sections in the previous section refer to the complete data set ineluding events from 
region III in whieh both D* and muon originate from the same b quark. Sinee this leads to a 
dilution of the correlation of quantities characterising the quark pair and the measured D*y pair, 
results for a smaller sample are given here, where both the heavy quark and the antiquark are 
tagged by either a D* or a muon (‘quark antiquark tag’). This is possible in eorrelation region 
IV (A® > 90° and Q{D*) ^ Q{y))- This region (see figure[T]) is dominated by cc pairs: the bb 
eontribution is 18% according to the two-dimensional fit. Due to migrations from eorrelation 
region III, approximately half of the bb eontribution is due to & — D*fi events in which the 
D*y pair comes from the same b quark (aeeording to the PYTHIA simulation). A visible cross 
section of a 2 q = 263 ± 48 ± 36 pb is measured^ after subtracting this fraction, while pb 

is expected from the FMNR calculations. 

In this data sample, the eorrelation of kinematie quantities reeonstrueted using the D* me¬ 
son and the muon to those of the quark pairs is good for Xg and y{D*y), while it is weaker for 
Pt{D*iS). Here Xg is the fraetion of the proton energy carried by the gluon in the hard inter¬ 
aetion, which is approximated by x°'^^ = {M{D*fx)Y/y s. The normalised differential cross 
section for x°Y^{D*iJi) is shown in figure |^. All QCD ealeulations (FMNR to FO and NFO, 
PYTHIA and CASCADE) give a reasonable description of the data. Figure 1^5 and e show 
the pT{D*y) and y{D*p) distributions of the D*p pair, respeetively, with the same model cal- 
eulations. The FO FMNR prediction for pt{D* p) is again too soft, as observed for the total 
sample (figure |5^), while the NFO FMNR prediction fits the data well. Although this sample 
should be sensitive to any transverse momentum of the ineoming gluon, the differenees between 
PYTHIA (eollinear faetorisation) and CASCADE {kt faetorisation) are small in the kinematie 
region studied. 


Conclusion 


A measurement of cc and bb photoproduction cross sections using the HI detector at HERA has 
been presented. For the majority of events both heavy quarks are tagged using a D* meson and a 

= 1/2 \n{E + Pz)/{E — P 2 ), where E and p^ are the energy and the z-component of the momentum of 
the D*y pair. 

^The index ‘2q’ is used for the cross sections in this section to distinguish them from those in the previous 
section. 
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muon as signatures. The separation of the eharm and beauty eontributions is possible due to the 
different eorrelations between the eharges and angles between the D* meson and the muon. The 
measured total eross seetion for eharm in the visible kinematie region is in agreement with the 
NLO QCD predietion, while the beauty eross seetion is higher than predieted. The ki nematic 
region of the latter is characterised by lower bb centre-of-mass energies than in most previous 
analyses, which require high momentum jets. Comparisons of the shapes of the measured 
differential distributions with QCD calculations including higher order effects show general 
agreement. Effects beyond the LO approximation are directly observed. In the kinematic region 
studied, effects due to kt factorisation, as implemented in CASCADE, are found to be small 
compared to the experimental errors. 
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I. Q{D*) = < 90° II. Q{D*) = Q(/i); > 90° 






Figure 2: a) Distribution of the mass difference AM = rnKniTs ~ f^K-wfor the total data sample. 
In b)-e) AM is shown in the four correlation regions, given by the relative charges of the D* 
and the muon and the azimuthal angle A® between them. The points represent the data (right 
charge combinations, RCC), the histogram indicates the observed wrong charge combinations 
(WCC) which are also used to fit the background. The solid lines in a) are the result of a one¬ 
dimensional fit, which gives the peak position p and the peak width a^M- The solid lines in 
b)-e) are results of a two-dimensional fit (see text). 
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Figure 3: Population of the four correlation regions I-IV obtained from the simultaneous likeli¬ 
hood fit in all correlation regions (two-dimensional fit, histogram). The resulting decomposition 
into charm and beauty contributions and the muon background is also shown. In all correlation 
regions the muon background is dominated by charm initiated events. The points with error 
bars are the results of one-dimensional fits of the AM distributions in each correlation region. 
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Figure 4: Normalised differential -D*/i cross sections as functions of the transverse momenta 
and the pseudorapidities of muons (a,c) and D* mesons (b,d). The data (points) are compared 
with the prediction of the NLO calculation FMNR and the LO-ME+PS QCD models PYTHIA 
and CASCADE. A beauty fraction of 45% as obtained from the measured cross sections is 
used in the calculations. The error bands for FMNR are obtained as described in tabled The 
PYTHIA h quark contribution is indicated separately. The experimental systematic uncertainties 
for the normalised distributions are negligible compared to the statistical errors. 
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Figure 5: Normalised dijferential D*^ cross sections for a,b) the transverse momenta 
c,d) the azimuthal angle difference A$, e) the invariant mass M{D*^) and f) the rapidity 
y{D*p) of the D*p-pairs. The data are compared to the prediction of the LO and NLO cal¬ 
culations FMNR (a,c,e,f) and to the Monte Carlo models PYTHIA and CASCADE (b,d,e,f). 
The error bands for FMNR are obtained as described in tabled A beauty fraction of 45% as 
obtained from the measured cross sections is used in the calculations. The FMNR (a,c) and 
PYTHIA (b,d,e,f) b quark contributions are indicated separately. The experimental systematic 
uncertainties for the normalised distributions are negligible compared to the statistical errors. 
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Figure 6: Normalised differential D*fi cross sections for a ‘quark antiquark tag’, charm domi¬ 
nated sample (approximately 10% bb quark contamination), where the D* and the p originate 
from different quarks. The data and predictions of the LO and NLO calculation FMNR and of 
the Monte Carlo generators PYTHIA and CASCADE are shown. The error bands for FMNR 
are obtained as described in tabled In c) the resolved excitation component of PYTHIA is in¬ 
dicated separately. The experimental systematic uncertainties for the normalised distributions 
are negligible compared to the statistical errors. 
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